Theory of the valley-valve effect in graphene nanoribbons 
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A potential step in a graphene nanoribbon with zigzag edges is shown to be an intrinsic source 
of intervalley scattering — no matter how smooth the step is on the scale of the lattice constant 
a. The valleys are coupled by a pair of localized states at the opposite edges, which act as an 
attractor/repellor for edge states propagating in valley K/K' . The relative displacement A along 
the ribbon of the localized states determines the conductance G. Our result G = (e^//!,)[l — cos(A'^7r + 
27rA/3a)] explains why the "valley- valve" effect (the blocking of the current by a p-n junction) 
depends on the parity of the number A'^ of carbon atoms across the ribbon. 

PACS numbers: 73.20.Fz, 73.23.-b, 73.40.Gk, 73.63.Nm 



I. INTRODUCTION 

The massless conduction electrons in a two- 
dimensional carbon lattice respond differently to 
an electric field than ordinary massive electrons do. 
Because the magnitude v of the velocity of a massless 
particle is independent of its energy, a massless electron 
moving along the field lines cannot be backscattered 
— since that would require f = at the turning 
point. The absence of backscattering was discovered in 
carbon nanotubes,— where it is responsible for the high 
conductivity in the presence of disorder. 

A graphene nanoribbon is essentially a carbon nano- 
tube that is cut open along the axis and flattened. One 
distinguishes armchair and zigzag nanotubes, depend- 
ing on whether the cut runs parallel or perpendicular to 
the carbon-carbon bonds. The edges of the nanoribbon 
fundamentally modify the ability of an electric field to 
backscatter electrons. As discovered in computer simu- 
lations by Wakabayashi and Aoki,^ a potential step in a 
zigzag nanoribbon blocks the current when it crosses the 
Fermi level, forming a p-n junction {— a junction of states 
in conduction and valence band). The current blocking 
was interpreted in Ref. ^ by analogy with the spin- valve 
effect in ferromagnetic junctions.* In this analogy the val- 
ley polarization in a zigzag nanoribbon plays the role of 
the spin polarization in a ferromagnet — hence the name 
"valley- valve" effect. 

It is the purpose of this paper to present a theory for 
this unusual phenomenon. A theory is urgently needed, 
because the analogy between spin valve and valley valve 
fails dramatically to explain the computer simulations of 
Fig. [T] The current blocking by the p-n junction turns 
out to depend on the parity of the number A'^ of atom 
rows in the ribbon. The current is blocked when N is 
even (zigzag configuration) , while it is not blocked when 
N is odd (anti-zigzag configuration, see Fig. [21). This 
even-odd difference (first noticed in connection with the 
quantum Hall effect^) is puzzling since zigzag and anti- 
zigzag nanoribbons are indistinguishable at the level of 
the Dirac equation,-'^ which is the wave equation that 
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FIG. 1: Conductance G of a zigzag nanoribbon contain- 
ing a potential step U = |(7o[l + tanh(22;/d)]. The red 
and blue curves are obtained by computer simulation of the 
tight-binding model of graphene, with parameters d = 10 a, 
Ef = 0.056 1, where a is the lattice constant and t is the 
nearest-neighbor hopping energy. Upon varying Uo the con- 
ductance switches abruptly to zero when the Fermi level Ef 
is crossed and a p-n junction is formed (red solid curve; the 
deviation from an ideally quantized step function is < 10~^). 
This "valley-valve" effect occurs only for an even number A'^ 
of carbon atom rows (zigzag configuration). When is odd 
(anti-zigzag configuration), the conductance remains fixed at 
2e^//i (blue dotted curve, again quantized within 10~^). 



governs the low-energy dynamics in graphene. 
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FIG. 2: Nanoribbons in the zigzag configuration (left panel, 
A'^ even) and in the anti-zigzag configuration (right panel, A'^ 
odd). In both cases the atoms at opposite edges belong to 
different sublattices (indicated by black and white dots). 



II. BREAKDOWN OF THE DIRAC EQUATION 
AT A POTENTIAL STEP 

The applicability of the Dirac equation rests on the 
assumption that a smooth potential step causes no in- 
tervalley scattering. As we now show, it is this assump- 
tion which fails in the p-n junction, breaking the analogy 
between spin valve and valley valve. In the spin valve, 
a spin-up electron incident on a ferromagnetic junction 
which only transmits spin-down is simply reflected as 
spin-up. The current blocking can therefore be under- 
stood without invoking spin-flip scattering. In the valley 
valve, however, an electron in valley K incident on a 
p-n junction which only transmits valley K' cannot be 
reflected in valley K. Both transmission and reflection 
require a switch of the valley from K to K' (see Fig. 
ini). We conclude that a p-n junction in a zigzag nano- 
ribbon is an intrinsic source of intervalley scattering. It 
does not matter how smooth the potential step might be, 
since the incoming and outgoing states are from different 
valleys, the scattering must switch valleys to preserve the 
current. 

As we have illustrated in Fig. [31 the source of inter- 
valley scattering is a pair of localized edge states at the 
p-n interface. It is well-known that the lowest mode in a 
zigzag nanoribbon is confined near the edges. ^ The trans- 
verse extension ^(e) ~ W/ In jeVF/?!-?;! of an edge state de- 
pends on the kinetic energy s = Ep — U. We define the 
p-n interface as the line where Ep--U{x, y) = 0. This line 
intersects the two edges at the points r± = {x±,y±), with 
2/+ = (|A'' — l)a/-\/3 = W and 2/_ = the y-coordinates 
of the row of atoms at the upper and lower edge, respec- 
tively. (Note that r± is well-defined also for a smooth 
interface.) Upon approaching the p-n interface, e de- 
creases from Ep to hence ^ decreases from £,{Ep) = 
to a minimal value of order of the lattice constant a. An 
electron incident on the p-n junction in valley K is there- 
fore attracted to a pair of localized edge states centered 
at r±. Their wave vector k spans the interval of order 
1/a between the valleys K and K' — thereby allowing 
for the intervalley scattering needed to repel the electron 
into valley K'. 




FIG. 3: Top panel; Zigzag nanoribbon containing a p-n in- 
terface from X- to x+ (dotted line). The spatial extension of 
edge states in the lowest mode is indicated by the grey areas. 
Incoming edge states are in valley K, while outgoing edge 
states are in valley K' . The arrows indicate the direction of 
propagation, in the conduction band (solid) and valence band 
(dashed). The corresponding dispersion relations are plotted 
in the lower panel. The localized (dispersionless) edge state, 
responsible for the intervalley scattering, is indicated in red. 



III. SCATTERING THEORY BEYOND THE 
DIRAC EQUATION 

Now that we have identified the mechanism for inter- 
valley scattering, we need to calculate the coupling of 
the propagating edge states to the pair of localized edge 
states, in order to determine whether an incident elec- 
tron is transmitted or reflected at the p-n interface. For 
this purpose we have developed a scattering theory based 
directly on the tight-binding Hamiltonian, 

Ho^t |n,m)(n',TO'|, (1) 

neighbors 

thereby going beyond the Dirac equation. The calcu- 
lation is outlined below, but we first present the result 
— which is remarkably simple: The transmission prob- 
ability T [and hence the conductance G = {2e'^ /h)T] is 
determined by the lateral displacement A — — X- of 
the localized states, according to 

T= i - icos(A^7r + 27rA/3a), (2) 

for ^ A. This is the central result of our paper. 

We have derived Eq. ([2]) by projecting the tight- 
binding Hamiltonian onto the pair of (nearly degenerate) 
lowest modes, and then solving a scattering problem in 
fc-space. As illustrated in Fig. [3] (lower panel), incoming 
and outgoing states have wave vectors near ki^ « An/Sa 
and fcout ~ 27r/3a, respectively. The unitary transfor- 
mation of an incoming state into an outgoing state is 
governed by the 2x2 transfer matrix M in the linear 
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relation 



= M{k,k')^{k'). 



(3) 



Here we have introduced the two-component wave func- 
tion \I>(fc) = (■0/l',V'fc ) fc-space. (For later use we also 
introduce the Pauli matrices at acting on the ± degree 
of freedom of the nearly degenerate lowest modes, with 
ctq the 2x2 unit matrix.) Once we know M, the scatter- 
ing matrix S = fiout^f (^out, Mn)f^in follows by a change 
of basis such that ^lx'^{kx) (with X labeling "in" or 
"out" ) has the first component in the conduction band 
(left end of the nanoribbon) and the second component 
in the valence band (right end of the nanoribbon). 

An analytical calculation is possible for W ^ £,o, when 
we can approximate the lowest modes ip^ = 2~^/^(-0^ ± 
V'f ) by^ 

^^(m,n) =C(fc)e™'='^/2[-2cos(fca/2)]"^'7r„+,„+i, (4) 
V^f (m,n) = C(fc)e™'='^/2[-2cos(fca/2)]'^""^n+™, (5) 



with C(fc) = \/— 1 — 2 cos ka a normalization factor. We 
have defined tt^ = 1 if p even and tt^ = if p odd. The 
integer n labels the row of atoms in the j/-direction and 
m labels the column of atoms in the a;-direction (see Fig. 
[2|). This approximation is accurate in the whole range 
(27r/3a, 47r/3a) of k, except within an interval of order 
1/W from the end points. The wave functions tp^, 
are edge states, extended either along the lower edge (on 
the A sublattice, indicated by black dots in Fig. [2]) or 
along the upper edge (on the B sublattice, white dots). 

The nearest-neigbor tight-binding Hamiltonian ([T]) is 
diagonal in the basis of the modes ■0^, with matrix ele- 
ments 

(fc, ±\Ho\k', ±) = ±e(k)a-^S{k - k'), (6) 
e{k) = 2tC{kf[-2cos{ka/2)]^. (7) 

Since e{7T/a — 6k) = {—l)^e(TT/a + 6k), the parity of N 
determines whether or not tpf, switches between conduc- 
tion and valence band as k crosses the point n/a. This 
band switch is at the origin of the parity dependence of 
the valley- valve effect, since it introduces a parity depen- 
dence of the matrices flin = <^o, ^out = crj^ that trans- 
form the transfer matrix into the scattering matrix. 
We model the p-n interface by a linear potential profile, 



Unm = + Uy{n~ N/2), 



(8) 



tilted by an angle 9 — arctan(| V3 ?7j^/[/a;). Upon pro- 
jection onto the two-component space spanned by ^'(fc), 
the Hamiltonian H = Hq + U becomes an integral kernel 
H{k, k') with a 2 X 2 matrix structure: 



H(k, k') = e{k)a-^6{k - k')a3 + iU^a''^—6(k - k')aQ 

ak 



The integral equation 



a I dk' H{k,k')^{k')^E^{k) (10) 

/27r/3a 

amounts to a system of two first order differential equa- 
tions: 

(^e(fc)a3 + lU^a'^ao-^ + ^NUy<Ji^ *(fc) = E-^{k). 

(11) 

This system gives directly an expression for the transfer 
matrix, 



M{k, k') = exp [i{k' - k)aE/U^]Texp 



e{q)a A 



I I dqfl{q) 
k' 



(12) 
(13) 



The scalar phase factor exp[i(fc' — k)aE/Ux\ has no ef- 
fect on the transmission probability, so we will omit it in 
what follows. The symbol T orders the operators in the 
exponent with respect to the variable q (from q = k at 
the left to q — k' at the right). The scattering matrix 
follows from 



5 = <M(fcout,fcin). 



(14) 



We may evaluate Eq. (|T2|) analytically ifW^A, be- 
cause then the integration interval can be separated into 
subintervals in which the contribution of one of the terms 
can be neglected. The calculation is described in App. 
El The result is 



M(fco 



exp[-i(7rA/3a)cri]e'" ' 



(15) 



+ ^Na~^Uy6{k-k')ai. 



(9) 



with a phase shift a = (— l)^a' that need not be deter- 
mined. Substitution into Eq. ^4]) yields the result ([2]) 
for the transmission probability T = 1 5*121^ ■ 

The regime A ^ W can be analyzed by a numerical 
evaluation of the integral HI]). The result, shown in Fig. 
[S] (solid curve), is that the conductance oscillations are 
damped for A > W. 



IV. COMPARISON WITH COMPUTER 
SIMULATIONS 

The current blocking (T = 0) obtained in the computer 
simulations of Refs.[3l3is the special case TV even, A = 0, 
corresponding to a zigzag configuration with potential U 
independent of y. In the anti-zigzag configuration {N 
odd) we have instead T 1, in accord with the simula- 
tions of Fig. [T] More generally, we can tilt the interface 
so that A 7^ 0. The simulations for a tilted p-n interface 
shown in Fig. [5] are well described by the analytical result 
(12), for A <C — 70 a. Note in particular the sum rule 
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FIG. 4: Conductance for a tilted p-n interface, with potential 
U = \Uo{l + tanli[2(x - ytan6l)/d]}, at fixed C/o = 2£f 
as a function of the relative displacement A = a;+ — 2;_ = 
W^tan^ of the intersection of the interface with the edges of 
the nanoribbon. The parameters are the same as in Fig. [T] 
which corresponds to A = 0. The data from this computer 
simulation is described by the analytical result ((2)). 




FIG. 5; Comparison between results of computer simulations 
(dots) and numerical evaluation of Eq. I|12p . The parameters 
are iV = 20 (solid line) and = 21 (dashed line), Uo = 
2Ef ^ 0.0058 1 and d = £f/;7^ = 100 a. 

G{N) + G{N + 1) w e^//i, first observed in the computer 
simulations of Ref. 5. 

For larger /S./W a phase shift appears and a reduc- 
tion of the amplitude of the oscillations, with G ~ for 
A^W . We compare the conductance calculated by nu- 
merical evaluation of the integral (fT2|) with the data from 
the computer simulations and find good agreement, see 
Fig. El 
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FIG. 6: The conductance of a tilted p-n interface with next- 
nearest neighbor hopping included {t' /t = 0.1). The parame- 
ters of the ribbon are Ef — 0.19 1, Uo = 0.16 1, and d = 100 a. 
The number of atoms across the ribbon is A'^ = 40 (solid line) 
and A'' = 41 (dashed line). 



W < AaT/Ep. The assumption can be relaxed in the 
case of a smooth p-n interface, because higher modes 
have an exponentially small transmission probability if 
the Fermi wavelength^ \p <^ d. 

Next-nearest-neighbor hopping was not included in the 
theory, and one might be concerned that it could mod- 
ify our result substantially because the edge states are 
then no longer dispersionlessjiS, We have found that this 
is actually not a relevant perturbation: Next-nearest- 
neigbor hopping (with hopping energy t') adds a term 
2t'(2 -t- coska)a~^5{k — k')(To to the projected Hamilto- 
nian ([9]). This is an irrelevant perturbation because its 
only effect is to multiply the transfer matrix by a 
scalar phase factor. 

As a check, we have repeated the computer simulations 
with the inclusion of next-nearest-neighbor hopping^ ^ in 
the tight-binding model (for the realistic ratio t' /t ~ 0.1). 
As shown in Fig. [SI the result ^ still applies for A <C W. 

Eq. ^ was derived for a linear potential profile U, but 
the derivation can be extended to include a smoothly 
varying potential landscape SU (smooth on the scale 
of the lattice constant). Electrostatic disorder there- 
fore affects the conductance only through the lateral dis- 
placement A of the points on the boundary at which 
U + SU - Ef =0. 

Edge disorder cannot be accounted for in this simple 
way, but in view of the small lateral extension of the lo- 
calized edge state we might not need a well-defined zigzag 
edge over long distances in order for Eq. ^ to apply. 



V. EXTENSIONS OF THE THEORY 

The theory presented so far can be extended in several 
ways. 

We have assumed that the width W of the nano- 
ribbon is sufficiently narrow that there is only a single 
propagating mode at the Fermi level, which requires^ 



VI. CONCLUSION 

In conclusion, we have presented a theory for the 
current blocking by a p-n junction in a zigzag nano- 
ribbon. The dependence on the parity of the number 
N of atoms across the ribbon, not noticed in earlier com- 
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puter simulations,—!^ is explained in terms of the parity 
of the lowest mode under a switch of sublattice: Incident 
and transmitted modes have opposite parity for N even, 
leading to complete reflection {G — 0), while they have 
the same parity for N odd, leading to complete trans- 
mission (G — 2e'^/h). A variation of the electrostatic 
potential in the direction transverse to the ribbon can 
invert the parity dependence of the conductance, while 
preserving the sum rule G{N) + G{N + 1) ~ 2e^/h. 

This switching behavior may have device applications, 
if the structure of the edges can be controlled (which 
is not the case in presently available samples). Even 
if such control is not forthcoming, the mechanism for 
current blocking proposed here can be operative in an 
uncontrolled way in disordered nanoribbons, producing 
highly resistive p-n interfaces at random positions along 
the ribbon. Conduction through the resulting series of 
weakly coupled regions would show an activated temper- 
ature dependence as a result of the Coulomb blockade;^ 
as observed experimentall y ^^'^'^ 
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APPENDIX A: EVALUATION OF THE 
TRANSFER MATRIX 

To evaluate the transfer matrix M(fcout,fcin) in the 
regime W ^ A we use the fact that the energy e{k) 
of the lowest modes decays exponentially ^ ex-p{—NaSk) 
away from the integration limits kin,kout [see Eq. (jT))]. 
We separate the integration in the momentum-ordered 
exponent l|12p into three intervals: 



M{kout, fcin) = Texp 



Texp 



■feout 

I I n{q)dq 



pk2 






i / Q{q)dq 


Texp 


i / n{q)dq 


Jki 







(Al) 



We choose A:in — A:i = fc2 — fcout ^ ^/W, such that |e(fci)| = 
|£(^2)| ^ Ux- Then the contribution of the term cTi A/2 in 
^l{q) to the integrals over the first and the third intervals 
is of order A/iy ^ 1, so that this term may be neglected. 
The contribution of the term <j^e{q)a/Ux to the integral 
over the second interval is of order a/W <C 1 so it can also 
be neglected. The three integrals can now be evaluated 
analytically, with the result: 



-M"(fcout, fcin) 



= e'""^ exp 



i{k2 



A 

^l)-^cri 



(A2) 



This is equivalent to Eq. (jisp since ^2 — fci = — 27r/3a 
0{l/W). 
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